


f 


ns Cain 






Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1991-03 


Optimal routing of military convoys through a 
road network 


Lee, Dong Keun 


Monterey, California. Naval Postgraduate School 
http://ndl.handle.net/10945/28430 
Copyright is reserved by the copyright owner 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
D U DLEY research materials and institutional publications created by the NPS community. 
get Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS'‘s first 
KNOX appointed — and published — scholarly author. 


LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


> o> om 

Ne = 
wee ore 
PR ee 


Lhe rs Pare en 
ee Ow nee 
2M Chey Gp a ra 

&. Sag 

ad . 


Wl 
i 


* 4 
sar f 34 at 


rr) 
te A 
= ae 
ee 
ee 

. 

ie 

Me 
a 


that 

















NAVAL POSTGRADUATE SCHOOL 
Monterey , California 








THESIS 


ae 


OPTIMAL ROUTING OF MILITARY CONVOYS 
THROUGH A ROAD NETWORK 


by 
Dong Keun Lee, 


March 1991 


Thesis Advisor R. Kevin Wood 





Approved for public release; distribution is unlimited. 


1253659 





lassified 


ity classtfication of this page 
REPORT DOCUMENTATION PAGE 


eport Security Classification Unclassified 1b Restrictive Markings 





ecurity Classification Authority 3 Distribution Availability of Report 

Yeclassification Downgrading Schedule Approved for public release; distnbution 1s unlimited. 
‘forming Organization Report Number(s) 5 Monitoring Organization Report Numberfs) 

vame of Performing Organization ob Office Symbol 7a Name of Monitoring Organization 

val Postgraduate School (if applicable) 55 Naval Postgraduate School 

.ddress (cin, State, and ZIP code) 7b Address (cid), state, and ZIP code) 

nterev, CA 93943-5000 Vlonterey. CA 93943-5000 





.ame of Funding Sponsoring Organization =| 8b Office Symbol! 9 Procurement instrument Identification Number 
(if applicable} 


.ddress (city, state, and ZIP code) 10 Source of Funding Numbers 


itle (include security classification) OPTIMAL ROUTING OF MILITARY CONVOYS THROUGH A ROAD NETWORK 
14 Date of Report (year, month, day) 


ersonal Author(s) Dong Keun Lee 
15 Page Count 
March 199] a5 


Type of Report {3b Time Covered 
ster’s Thesis From To 

upplementary Notation [he views expressed in this thesis are those of the author and do not reflect the official policy or po- 
yn of the Department of Defense or the U.S. Government. 

rosauu Codes 18 Subject Terms (continue on reverse if necessary and Identify by block number) 


Multycommodity, Dynamic, [Transportation network. 












.Ostract (continue on reverse if necessary and identify by block number) 

To wage a successful campaign, military units and materiel must be in position by the designated tume. This thesis models 
problem of moving military units and matenel in convoys through a road network as mathematical programming models. 
articular, two models, linear and integer, are investigated. Both models belong to the class of multicommodity, dynamic 
isshipment network problems. Based on prototypic GAMS implementations, they provide essentially the same answer. 
wever, the linear model ts easier to construct, takes less tume to solve and allows for more flexible convoy routing. 


Distribution Availabilty of Abstract 21 Abstract Security Classification 
unclassified unlimited CJ same as report CL} DTIC users Unclassified 
Name of Responsible Individual 226 Telephone (include Area code) 22c Office Symbol 
Kevin Wood 08) 389 


FORMI 1473,84 MAR 83 APR edition may be used until exhausted security classification of this page 
All other editions are obsolete Oe 


Unclassitied 


pao 


Approved for public release; distribution 1s unlimited. 


Optimal Routing of Military Convoys 
Through a Road Network 


by 
Dong Keun Lee 
Major, Republic of Korea Army 


B.S., Korea Military Academy, 1982 


Submitted in partial fulfillment of the 
re “lirements for the decree of 


MASTER OF SCIENCE IN OPERATION RESEARCH 
from the 


NAVAL POSTGRADUATE SCHOOL 
March 1991 


ABSTRACT 


To wage a successful campaign, military units and materiel must be in position by 
the designated time. This thesis models the problem of moving mulitarv units and 
materiel in convoys through a road network as mathematical programming models. In 
particular. two models, linear and integer, are investigated. Both models belong to the 
class of multicommodity, dynamic transshipment network problems. Based on 
prototypic GAMS implementations, they provide essentially the same answer. However. 
the linear model ts easier to construct, takes less time to solve and allows for more flex- 


ible convoy routing. 
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I. INTRODUCTION 


During the initial phase of a war, a large number of troops and materiel must Si- 
multaneously move up to the front line in convoys. Routing of these convoys to their 
final destinations is inherently difficult because many possibilities exist. When routing 
is done nonoptimally, certain routes may become clogged while others are left unused 
and this may lead to delays detrimental to military objectives. Therefore, it is the goal 
of this thesis to develop and solve the problem of routing convoys to their destinations 
by designated times. 

1. AIRLAND BATTLE 

The concept of the Airland Battle in current military doctrine is a method to 
defeat a large armored force which attacks by echelon, through a narrow breach sector 
in am opponent's front. The Airland Battle requires complex integration of fire support 
weapons, Army and Air Force aviation assets and ground maneuver forces. Doctrine 
calls for the enemy to be defeated by first stopping his forward elements, second, by 
defeating his reserve echelons before their combat value can influence the battle, and 
third, by attacking his resources to prevent reconstitution of previously defeated forces. 

Defeating the armored force echelons at the proposed breach point will require 
more troops and war niateriel than in past conflicts. Most of the troops and materiel 
will arrive at the front by land convoys. Thus, the problem of effectively moving troops 
and war materiel through the available road network in a short period of tme 1s even 
more important than before. 

2. PROBLEM SCOPE AND GOAL 

In this thesis, the term “unit” refers to a group of troops and materiel which 
will move through a road network as a single, coherent convoy. A convov consists of 
at least 6 vehicles moving at the same time or 10 vehicles moving within a l-hour period 
under a single commander, over the same route and in the same direction [Ref. I : p.5-1]. 
In practice, there can be several hundred vehicles in a single convoy and each of several 
convoys must arrive at its destination in a limited amount of time. 

This thesis describes two mathematical programming models which determine 
routes for multiple convoys to arrive at their destination by the designated times. To 
determine the effectiveness of the models, they are implemented in GAMS (General Al- 


gebraic Modeling Systein) and their solutions are compared. 


3. THE MODELS AND SOLUTION PROCEDURES 

To simplify the presentation, a concrete scenario is used to demonstrate the 
models throughout this thesis. In the scenario, there are three units at three separate 
staging areas behind the battle front. Each unit must move through the road network 
to its own individual supply point to take on supplies and then move up to a common 
front-line location. A certain nunimum amount of delay at supply points is incurred 
since some time 1s required to load supplies. The two models investigated in this thesis 
have the same underlving network structure in that the units movements are represented 
as flows through a road network. However, they have different assumptions on convoy 
moveinents. 

The first model is an integer program which assumes that a unit consists of only 
one convoy and all convoys are of the same column length which is the length from the 
head of the convov to the end of convov. Moreover, all arcs in the network have the 
same capacity which allows only one convoy on any road link at any time. This sim- 
plifies the coordination of troop movements. 

The second model is a linear program which allows units to move in convoys 
with different lengths. Each unit can be split up into smaller (sub)convoys, each taking 
different routes. Also, arcs can now have different capacities and can hold several 
(sub)convoys from different units. Although this model potentially uses the road net- 
work more efficiently, it is more complicated to administer and coordinate. This model 
can also be viewed as a multicommodity extension of the building evacuation models in 
which occupants of a large building must be evacuated in a short amount of time. [Refs. 
p53) 

It is important to note that in both models units starting from a different staging 
area must be treated as separate commodities in a network flow formulation. Otherwise, 
it cannot be guaranteed that a complete unit will pass through its own supply point and 
end up at the correct destinations [Refs. 2,3,4,5]. In addition, the network must be dy- 
namic in that the basic road network is expanded over time to represent the fact that a 
certain amount of time is necessary to transit various links in the network. 

4. OUTLINE 

Chapter II describes the concept of a multicommodity dvnamic network model. 
Then, the two basic models are presented in Chapter III. Finally, the computational 
results, conclusions and recommendation for further research are given in Chapter IV. 


For completeness, the GAMS programs are listed in the appendices. 


1. MULTICOMMODITY DYNAMIC ROAD NETWORK 


This chapter describes how the basic structure of a road network is transformed into 
a multicommodity dynamic network for purposes of modeling convoy routing problems. 
[In the basic road network, nodes represent road junctions, troop staging areas, supplv 
points and the front line. Arcs represent segments of the road network. In the mulu- 
conimoditv dynanuc network, nodes, for the most part, correspond to nodes in the ori- 
ginal network but are replicated over multiple time periods. Arcs in the dynamic 
network correspond to movements of units through time and space. 
1 DIRECTED AND UNDIRECTED GRAPHS 
Let G =(V,E) be a graph which represents a road network. V is the set of nodes 
and £ is the set of arcs where each arc ee E consists of a pair (u,v) where ue V and 
ve V. A graph may either be directed or undirected [Ref. 6 : p.198, : p.230]. In an 
undirected graph any arc (u,v) is taken as an unordered pair while in a directed graph the 
pair is taken as ordered. In a directed graph an arc (u,v) is often expressed by u > v and 
u is called the rai of the arc and v is called the head of the arc. A road network can be 
represented bv a directed graph with one-way road segments being represented by single 
arcs and two-way segments represented by two arcs in anti-parallel. While most roads 
networks consist mainly of two-wav segments, for the purposes of moving units up to a 
front-line, most of the segments can be considered to be one-way, oriented in the direc- 
tion of the front. 
2. SINGLE AND MULTICOMMODITY DYNAMIC NETWORKS 
Consider now a homogeneous, infinitely divisible commodity such as 
“materiel” measured in units of, say, tons moving through a road network G=(V,E) . 
G is called the static network since its structure does not include a reference to time. 
However, since each arc ee E in this network will require some finite time to transit It, 
an integer transit time t(e) is associated with the arc. 
Consider the network G;=(V,,E,), the T-time expanded network obtained from 
G =(V, E) as follows: 


Ver= {uu € V, = Pare Be 


Here u, is the t-th time copy of node we V. Simularly, the arc set E; is given 


E-= (jy, jie = (0) es — tae et = ee le 
WG, ie Vor 1,2... 7 aa 


The network G, 1s called the dynamic netvork associated with G . The arcs 
(u,, U,.,)€ E; are called the holdover arcs. Traversing such an arc represents materiel 
pausing at node uw from time period s to time periodr+ 1. 

The arcs (u,, v,) are called the movement arcs. They represent the movement of 
materiel from one node in the road network to another. This movemient starts at time # 
at node uw, and terminates at node v at time s . Associated with each holdover arc 
(t,, U,.) , 1S A Capacitv which represents, for example, the tons of materiel which can 
pause at node uw in the road network. Associated with each movement arc (u,, ¥,) , Is a 
Capacity which represents the tons of materiel which can be moved from uw to v in 
(s — t) time periods. 

Using standard flow balance constraints associated with the dynamic network 
[Ref. 3 : pp.98-99], and defining supplies and demands, it 1s then possible to model the 
flow of materiel through time and space through the road network. This is not sufficient 
for the purposes of this thesis, however, since multiple commodities, 1.e., different mill- 
tary units, must be distinguished. Consequently, a multicommodity variant of the dy- 
namic network model must be formulated. 

In the multicommodity dynamic network, the basic dynamic network 1s repli- 
cated for each commodity. Each commodity network has its own supplies and demand. 
Then, as in standard multicommodity network flow problems [Ref. 5], joint capacity 
constraints are placed across the commoc. ies for each set of analog-:us arcs. This then 
describes the setup necessary for the two models described in the next chapter. Some 
restrictions will have to be added however, since, in the first model, it will not be as- 
sumed that the commodities, i.e., military units, are infinitely divisible. Furthermore, 
certain modifications to the dynamic networks will be necessary to accommodate differ- 
ent objective functions and constraints. 

3. TIME-EXPANDED NETWORKS 

The following discussion describes some of the issues associated with practical 
use of the dynamic networks. Consider a road which is composed of an arc and two 
nodes. Figure 1 shows the static and dynamic network representation of this road net- 


work. 


STATIC DYNAMIC 





Figure 1. Static and Dynamic Networks 


I-or each static directed arc, say froin static node wu to static node v , and having 
traversal time | and for every integer ¢ between 0 and T—1 construct a (directed) 
movement arc in the dynamic network from copy ¢ of node u to copv 1+ 1 of node r. 
It is easy to verify that if the static model has » nodes and a arcs, and the dynamic model 
has 7 time periods, then (7+ a)T is an upper bound on the number of arcs tn the dy- 
naimc model, and 2(7'+ 1) is an upper bound on the number of nodes of the dvnaimic 
model. These upper bounds can often be decreased substantially by deleting “Inessen- 
tial” arcs and nodes in the dynamic model, 1.e., arcs and nodes not lving in at least one 
directed path from copy 0 of some tail node to copy 7 of some head node. 

Let R be the maximum transportation planning time, and [et t be the length 
of a time period. Then, 7, the number of time periods in the model is given by 7 = K/r 
assuming exact divisibility. Thus, the dynamic model can have as many as n(R/t + 1) 
nodes, in which case its computational tractability will be inversely proportional to the 
magnitude of ft. 

Ideally, supposing all arc traversal times to be integers originally, a reasonable 
choice of t is the greatest common divisor of all the traversal times. Unfortunately the 
greatest common divisor may be onc, in which case it may well be necessary to alter 


some of the traversal times in order to find an acceptably large greatest commion divisor. 


In the view of the emergency nature of the convoy routing problem one mav wish to 
make all such alterations larger, in order to be assured that the model will not underes- 


timate the minimum transportation time. [Ref. 2: pp.90-91.] 


It. MODEL DEVELOPMENT AND FORMULATION 


Time-minimizing transportation routes and schedules are concerned with determin- 
ing the nunimum time required for transporting a set of goods from given supply points 
to given demand points. The importance of such problems arises when certain con- 
sumers urgently require a set of commodities to be supplied from fixed distribution cen- 
ters. The convoy routing problems are examples or at least modified examples of such 
time-minimizing problems. This chapter details two models for convoy routing prob- 
lems. 

I. SCENARIO 

Consider a country under defense of a border point. The defense must be met 
by a given number of units and their equipment. It will be assumed that there are three 
units located at three separate points along with their equipment but not all of their 
supplies. Each unit must first move forward through the road network to its own unique 
supply point where it takes on supplies. A certain amount of time is necessary to take 
on supplies; if less time is spent a unit would have to move to the front line without its 
full complement of supplies making the unit less than 100% effective. After taking on 
supplies, the units must move up to a single front line location to arrive as early as 
possible or at least bv a specified time. For this transportation operation to be suc- 
cessful, the sum of the time to move, the time to take on the equipment and unknown 
delay time should be less than the allowed operation time. 

In these models, the arrival time at the destination and the sojourn time at the 
supply points are critical. If the units do not spend enough time at their respective 
supply points thev will not have sufficient time to take on supplies. If the units do not 
move to the destination within the required time, a coordinated attack at the breach 
point cannot be launched and the operation may fail. Some conflict in these objectives 
may arise. 

Consequently, two different objective types for Model | will be considered. In 
the one type, it is assumed that the units should arrive at the destination as early as 
possible so as to have as much time as possible at the front line to prepare for battle. 
However, they must spend a minimum amount of time at the supply points. [n the other 
type, it is assumed that the units need only arrive at the destination by a specified time 


and they should spend as much time as possible at the supply points taking on supplies 


and perhaps making other preparations. It will be assumed that the length of a convoy, 
called its “columr. [ength” ,is fixed. 

Two oifferent column lengths for Model 2 will be considered with one objective, 
to minimize the arrival time at the destination. First, the number of vehicles in a convoy 
is equal to arc Capacity in order to compare the solution to Model 1. Second, the 
number of vehicles in a convoy 1s made greater than arc capacity so that the convov 
must be broken down into at least two small subconvoyvs. In addition, the capacity of 
one arc 1s reduced to 0 to represent a situation in which part of the road network is de- 
stroved bv enemy attack. 

2, FIXED COLUMN LENGTH, MODEL I 

In the first model, Model 1, it 1s assumed that the column length of each unit 
corresponds to a fixed time period t. Thus, if a particular arc in the road network re- 
quires one time period to traverse and the column length 1s two time periods. it will take 
two time periods for the unit to clear the arc. [t is assumed that the joint capacity of 
anv arc and node, except the staging nodes, supply point nodes and the destination 
nodes are 1. Decision variables in this model will be binary indicating whether or not a 
unit 1$ transiting a particular are or sojourning at a particular node. Thus, this model 
iS an integer programming model (IP). 

a. Procedure 

Figure 2 shows the road network which 1s used to test Model 1. The num- 
ber on each arc is its traversal time. To successfully model fixed column lengths it 1s 
necessary to split the arcs in the network into multiple arcs in series so that each has a 
length equal to a fixed fraction of the column length: I, 1/Z,1/3, ete. Thus, if the lengun 
of each arc were 1/3 of the column length it would take 3 fime periods for the column 
to clear the arc, Clearlv, the shorter the length of each arc the better for accuracy but 
this, of course, leads to larger problems. For simplicity and computational tractability, 
it is assumed that each arc can be split into lengths which are equal to the column length 
(The column length 1s fairly short compared to actual arc lengths.). Figure 3 shows the 
modified static network in which each arc has the same length. . 

Once the static network is obtained it can be expanded into a dynamic net- 
work over a specified number of time periods. The number of time periods should be 
as few as possible to lead to a model which is as small as possible. Figure 4 displays the 
dynamic network. 

It may be necessary for one or more units to wait to allow another unit to 


pass. Thus, it is sensible to allow holdover arcs on all nodes. However, this leads to a 





source supply 
unit-f f 6 
unll-2 2 1! 
unt{-3 3 13 


Figure 2. Road Network 
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Figure 3. Modified Network 


verv large number of arcs and to avoid this it is assuined that any waiting will be done 
at the source nodes, the supply point nodes or at the destination node. Thus, only 7 
nodes in the road network need holdover arcs. 

The number of nodes and arcs in the dvnamic network can be reduced fur- 
ther. Every node 1s expanded by time period. Ilowever, there can be no flow in the first 
time period on any arcs except those adjacent to the source nodes. Thus, many arcs can 
be deleted. figure 5 shows the reduced dynamic network. It would be desirable to 
eliminate all arcs which can never have any flow on them because they are correspond 
to an early time period too far away from any source node. This shows on the night 
upper right of Figure 5. Also, there can be no flow on arcs too far away from the des- 
tination at later time periods. Thus, additional arcs could be deleted as shown in the 
lower left portion of Figure 5. The destination node does not need the holdover arc 
because its capacity is unlimited. The next step is to formulate this problem. The for- 
mulation assumes that the column length equals the arc lengths. The modifications 
necessarv to handle longer columns are discussed in Chapter IV. 

b. MNlodel 1 Formulation 


1. Indices 
® u,v=1,2,...,N nodes 
ean — [.2,...,K units 


Seg — li 2 me periods 


2. Data 
a, = 1 if there exists an arc between node u 
and node v at time period f 
0 otherwise 
e mintime,, minimum time to take on 
equipment at supply point uw for unit kK. 
@ = arccapacity,,, capacity of arc between u and v 


¢ nodecapacity,, capacity of node w at time period ¢ 
3. Decision Variables 


© Xun = 1 if unit & traverses arc (u,v) 
at time ft 
QO otherwise 


7 total amount of time to transport 
4. Formulation 


e The objective function ts 
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Modified Dynamic Network for Model | 


Figure 5. 


min »y » » L Aynr—y Xeuvji-1 Where N is the destination node (3.130) 


kK {ou 


e Subject to 


> Qyyr Xenyur= 1 for Vk, and v the source node for k,t=1 (3182) 
ld 


s » Quye Xpuyr—1 + » Aur Xk = 9 for Vk, t, and v (3.1.3) 
u td 


y » Quy Xeyvp=—l for WV k and where N is the destination node. (3.1.4) 
fe u 


»y Qivr Xiqjyp SS QNCCADGEIY,.,, [OO ama al (3.13) 
k 


» » Quut Xkuyr 2 Mintime,, for Vk, and u the supply point for k (3.1.6; 
uot 


» » Gy; Xiuy: = MCGCCAPaciy a or N veer (3.0 
k u 


The GAMS code for Model 1 1s given in appendix A. 

The decision variable x,,,, establish a flow amount from node u to v at time 
period ¢ for unit &. If the objective is to minimize the average arrival time of the units 
at the destination, the objective function is 3.1.1. However, if the objective is to maxi- 


mize the average time spent at the supply points , the objective function 1s 


max y » » Quut Xkuur Where U is the set of the supply points (3.53) 
kK ue fF 


Constraints 3.1.2,3.1.3, and 3.1.4 are the flow balance equations. Con- 
Straints 3.1.2 specify one unit of “supply” at the source node for each unit &. Constraints 
3.1.3 are balance equations for intermediate nodes,1.e. nodes other than the source and 
destination nodes. Constraints 3.1.4 for node N, the destination node, state that each 
unit must arrive at the destination in some allowable time period. Constraints 3.1.5 are 
joint capacities for each arc and time period. For the implementation of this thesis all 
Capacities are 1. Constraints 3.1.6 enforce a minimum sojourn time at the supplv point 


for each unit. Constraints 3.1.7 are the joint node capacity constraints corresponding to 


capacitated holdover arcs. The destination node has infinite capacity and 1s not included 
Nere: 
3. VARIABLE COLUMN LENGTH : MODEL 2 

In Model 2 it is assumed that any convoy may be stretched out and intermingled 
with other convoys subject to capacity limitations such as the maximum number of ve- 
hicles that can pass through an arc at any one time. However, the identity of a unit 
must be maintained since it has its own supply point which must be visited. In this 
model, the decision variables, the number of vehicles traversing an are at any one time, 
are allowed to be continuous. Thus this model is a linear program (LP). 

a. Procedure 

Figure 2 shows the static road network with arc traversal times. In Model 

1, the static network is first modified so that each arc has the same length as a column. 
Model 2 does not enforce a fixed column length so that the static network can be directly 
expanded into a dynamic network with some arcs spanning | unit of time, some 2 units 
of time, some 3 units of time and so on. This dynamic network is shown in figure 7. 

In model 2 , any unit may be broken down into several smaller units which 
mav use different routes. If the same constraints are used as in Model | to force a 
sojourn at the supply points some part of a unit may skip the supply point entirely and 
another part may spend extra time at supply point to satisfy the constraint. To avoid 
this, holdover arcs at the supplv points are used which have the same length as the time 
needed to take on equipment. This 1s show 1s shown 1n figure 7. 

In figure 8, nodes and arcs which can have no flow are eliminated. In ad- 
dition, junction nodes have holdover arcs since some delay may occur there; such arcs 
have restricted capacity, too. The destination node doesn’t need holdover arcs because 
its capacity 1s unlimited. 

b. Model 2 Formulation 
1. Indices 
e u,v=1,2,...,.N nodes 
e k= 1,2,...,K units 


© t,s=1,2,...,/ time periods 


2. Data 
® Qa,= | if there exists an arc between node u and 
node v from time period f to time period s 
QO otherwise 
¢ vehicles, vehicles in unit k 


5 


OO 0 O 
OOO0 O 
CRCRORORe 
Oe C6) C7 @ 


© SE" O2@ 7 © 
YQNaAaN A 


Dynamic Network 


Figure 7. 


ee | maximum time for the operation 
® arccapacity,,,, Capacity of arc (u,v) between periods s and ¢ 
°® nodecapacity,, capacitv of node w in time period r 

3. Decision Variables 


© Xnes= number of vehicles for unit & traverses arc (u,v) 
from time period f to time period s 


9ie7 total amount of time to transport 
4. Formulation 


e The objective function is 


min » » y Ss S Anis Xkunrs Where N is the destination node. (3.2515 
Piet ou Gas 


e Subject to 


>> Avurs Xkyyrs = Vehicles, for V k, v the source node fork andr=1 = (3.2.2) 
5 u 


= y yee Spy » » Ayre Xtuts =O for Vk, v,t (3.2.3) 
u 5 u 5 


» yy Dy ane XpuNse = — vehicles, for Yk, N the destination node (3.2.4) 
u I g 


» Giver Xtyver ear c CCD AGI as OF MVE, «S. 1 (3.2. 
k 

ie > Ayy,st Xkuv st = ever) Xkv,v, eer) for k,t (3.2.6) 
u 5 

uv 


where r 1s sojourn time and y, is supply point v for unit & 


» » » Ayysr Xkuvst S NOdecapacily,, for Vv,t (3:22) 
k u 5 
» y, » Qyvis Xkuves = vehicles, for V k, v the supply point for k (3.2.5) 
u I S 

uxV 


The objective function 3.2.1 minimizes the average arrival time to the des- 
tination node over all units. There ts no variant of this model, as in Model 1, in which 


the objective is to maximize the sojourn times at the supply points. Constraints 3.2.2 


are the balance equations in the first time period at the source nodes. Constraints 3.2.3 
are balance equations for intermediate nodes and constraints 3.2.4 are for the destination 
node for each unit. Constraints 3.2.5 are the joint arc capacities for arcs (u,v) between 
time periods s and ¢. Constraints 3.2.6 enforce the necessary sojourn time at the supply 
point v for each unit. Constraints 3.2.7 enforce node capacities. Constraints 3.2.8 re- 


quires that for each unit all its vehicles pass through the unit’s supply point. 
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IV. IMPLEMENTATION, RESULTS AND CONCLUSIONS 


This chapter describes the results from the GAMS implementation of Models | and 
2 on an IBM 3033AP at the Naval Postgraduate School(NPS). The GAMS complier 


at NPS uses BDMLP and ZOOM to solve linear and integer programs, respectively. 


Il MODEL I: FIXED COLUMN LENGTH 

Two versions of this model, Models 1.1 and 1.2, were implemented. Model 1.1 
nunimizes the average arrival time of the units at the destination node and Model 1.2 
maximizes the sojourn time at the supply points. 

The road network contains 44 nodes and bv solving three shortest path prob- 
lems, one for each of the 3 units, it is determined that at least 24 time periods are needed 
for feasibilitv. In fact, 25 time periods were used and this proved sufficient. If all possi- 
ble nodes and arcs are generated, the complete integer program would contain 145,200 
decision variables (arcs) and 52,803 constraints. However, some of the variables (arcs) 
and flow balance constraints may discarded since they do not affect the optimal solution. 


The model statistics for models 1.1 and 1.2 are summarized in Table 1. 





2867 2025 
7912 0-1 variables 2022 


ables : 


Model 1.1 21.30 sec 21.40 sec 
Model 1.2 21.35 sec 


Wales: MODEL STATISTICS FOR MODELS 1.1 AND 1.2 









The optimal paths for Models 1.1 and 1.2 are shown in Figure 8 and 9, respec- 
tively. Model 1.1 has no delays since unit I arrives in 13 time periods, unit 2 in 15 time 
periods and unit 3 in 24 time periods all of which are the times which can be achieved 
if no joint capacity constraints are enforced. In the solution to Model 1.2, the units 
delav at the supply points as possible as and all arrive by time period 24 as required. The 


route taken by unit 2 is different than that taken in Model 1.1. 
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2, MODEL2: VARIABLE COLUMN LENGTH 
Simular to the fixed convoy length model, two versions of this model, Models 
2.1 and 2.2, were implemented. To compare Models | and 2, the capacities for all nodes 
and arcs in Model 2.1 are set to the length of the convoy, which is assumed to be 30., 
The capacities of nodes and arcs in Model 2.2 are the same as in Model 2.1; however, 
the convoy length is now 36. The road network for Model 2 is the same as the one in 


Model 1. The statistics for Models 2.1 and 2.2 are summarized in Table 2. 


Vode 2 This see 
Model 23 1.07 sea 1.16 seq 


dable2. MODEL STATISTICS FOR MODELS 24 AND 2:2 









Figure 10 and 11 display the optimal solutions of Models 2.1 and 2.2, respec- 
tively. In the solution of Model 2.1, there is no delay experienced bv any of the units 
because all arc and node capacities are sufficiently large to handle all the traffic. In the 
solution of Model 2.2, however, there is some delay. In particular, a fraction of unit 2 
(6 vehicles) uses the arc between nodes 5 and 6 in period 3 while a fraction of unit | (24 
vehicles) uses same arc at same time period. So, a fraction of unit 1 (24 vehicles) goes 
to node 5 and the other (6 vehicles) stays at node 4 for | time period. Unit I is then 
broken down into two components one of which contains 24 vehicles while the other 
contains 12 vehicles. The two components of unit | still use the same route to the des- 
tination node; part of the unit just lags behind the other. In fact, each unit uses a single 
route. 

If some arc capacities are reduced, there might be more delays and the units 
might split into components following different routes. In fact, if the capacity of arc (9, 
13) between time periods 14 and 16 is reduced to 0 then unit | splits up at node 6 and 
part of the unit reaches the destination via nodes 9,13,15, and 16 while the other part 
uses nodes 8,14,15, and 16; the other units use the same routes as in Figure I 1. 

3. COMPARISION OF BOTH MODELS 
Below, we list the advantages and disadvantages of the two models. 
1. Model 1: Advantages 


a. The objective function is straightforward. 


2k 





time-period 


Ch 


(0) __ untt-2 


source supply sink 





unit 1 1 7 44 

| 
unit 2 2 23 44 
Q 2 unit 3 3 27 44 


Figure 8. Optimal Paths For Model 1.1 


22 





time-pertod 


ee om =_— unit- 


———$p— IE -2 


ak oe ee NIE =D 
| Ca: 


source supply sink 
unit f / 7 44 

@ unlit 2 2 23 44 
@ a @ unit J 2 | ra 44 


Figure 9. Optimal Paths For Model !.2 
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Figure 10. 


Optimal Paths For Model 2.1 
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Figure 11. Optimal Paths For Model 2.2 
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b. It is easy to modify the minimum time required to take on supplies at the supplv 
points: Just change the data element runtime, . 


2. Model I: Disadvantages 


a. The length of all arcs should be the same and must be related to the column 
length. This decreases the accuracy of model if there is a large difference be- 
tween columin length and the length of the shortest arc. 


b. The number of nodes and arcs in the static network is increased by the necessity 
to modify the length of all arcs to be the same. 


c. It is hard to change the column length. For instance, suppose that the length 
of each arc remains one time unit but that the column length is two time units. 
Then, x;,,,, becomes x,,,,,, [hat 1s, we must consider that a unit is covering two 
arcs (u,v) and (v.w) at any one time. This leads to many more variables in the 
formulation. Furthermore, the balance constraints become much more compli- 
cated and numerous. 


d. Because Model | 1s an integer program, not an LP as in Model 2, computational 
times are long. 


e. Model 1 has more variables and equations than Model 2 after applying re- 
ductions. 


3. Model 2 : Advantages 


a. The dynamic network for this model is easier to construct than the analogous 
network for Model | since it is not necessary to modif¥ the initial static network 
so that all arc lengths are equal to the column length. 


b. This model has more flexibility in that it is easv to change the column length 
Without modification of network structure. (Just change the supplies and de- 
mands.) 


c. Model 2 has fewer variables and equations than Model | after reductions. 


d. Model 2 is an LP and computational time is much shorter t in for Model 1 
which is an IP. : 


4. Model 2: Disadvantages 


a. It is hard to modify the time to take on equipment since it is given by network 
structure. 


b. There does not appear to be a variant of Model 2 analogous to Model 1.2 in 
which the average time spent at the supplv points 1s maximized. 
4. CONCLUSIONS AND RECOMMENDATION 
Based on the GAMS implementation using the example scenario, Model 2 pro- 
vides essentially the same answer as Model | using considerably less cpu time. In addi- 
tion, Model 2 1s more flexible and can handle larger problems. )n the other hand, 


ta 
Model 2 has one disadvantage in that it can not maximize the ~ spent at supply 
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points. Flowever, one could use Model 2 to find the maximum time by varving the 
length of sojourn at supply points. 

Future work in this area should include applying Model 2 to larger scenarios 
and extending the model to the case which includes the efrects of congestion. This has 
been done, for instance. in building evacuation models where the flow of people through 
an exit route can be reduced when the densitv of people becomes too high. This non- 
linear effect can be included, at least approximately, in a modified linear programming 
model [Ref. 3]. 
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APPENDIX A. GAMS PROGRAM FOR MODEL 1 


STITLE MIXED INTEGER PROGRAM FOR A TIME PERIOD COLUMN LENGTH 

“THESIS eMODEL 

*MAJ DONG KEUN,) LEE 

*MODEL: MULTICOMMODITY DYNAMIC TRANSPORTATION NETWORK TO MINIMIZE THE 
ue TRANSPORTING TIME. 


SETS 
N NODES /N1*N4&4/ 
T TIME EERIODS e022) 
K #7 OF TROOPS ~~) Kick ae, 
ALIAS(N,M); 


PARAMETERS 
A(N,M,T) ROAD NETWORKBETWEEN N AND M AT TIME T FOR TROOP K 

/N1. (N1,N4).T1*T11 
N2. (N2,N4,N5). T1*T11 
N3. (N3,N6). T1*T11 
N4.N5. T2*T12 
NS. (N7,N8). T2*T12 
N6. N8. T2*T12 
N7. (N7,N11,N9). T3*T15 
N8.(N13. +). T3*T16 
N9.(N15,N10). T3*T16 
N10. N16. T4*T18 
N11.N12. T4*T16 
N12.N14. T5*T16 
N13.N18. T4*T14 
N14. N18. T6*T18 
N15.N17. T5*T18 
N16.N19. T6*T18 
N17.N21. T6*T20 
N18. N27. T5*T15 
N19. (N20,N26). T7*T20 
N20. N23. T8*T20 
N21. N22. T7*T20 
N22. N28. T8*T20 
N23. (N24 ,N30,N23). T9*T19 
N24. N25. T10*T20 
N25.N26. T11*T21 
N26. N28. T12*T22 
N27. (N29 ,N27). T6*T20 
N28. N29. T9*T20 
N29.N31. T7*T23 
N30. N32. T10*T20 
N31. N34. T8*T23 
N32. N33. T11*T21 
N33. N35. T12*T22 
N34. N36. T9*T24 
N35. N37. T13*T23 
N36. N39. T14*T24 
N37. N38. T13*T23 


See ee ee ee eee eee ee eee eee en ee ee oe ee de ee el oe 
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N38. N40. T13*T23 
N39. N4&4. T1i5*T24 
N40.N41. T12*T22 
N41.N42. T13*T23 
N42.N43. T14*T24 
N43. N44. T15*T24 


el oe ee ee ee on 


i 


MP(M) SOURCE SUPPLY POINT AND SINK NODE 
/N1 1 
N2 1 
N31 
N7 1 
N23 1 
N27 1 
N44 1 


MIDICK,N,N) SUPPLY POINT NFOR TROOP K 


/K1.N27.N27 1 
K2.N23.N23 1 
K3.N7.N7 1 


S(K,N,T) SOURCE AND SINK NODE FOR UNIT K AT TIME T 
/K1.N3.T1 1 
RZ N2e Ti 1 
K3.N1.T1 1 /; 


SCALAR MINTIM TIME TO TAKE ON EQUIPMENT AT SUPPLY POINT /5/; 
SCALAR ARCCAP ARC CAPACITY /1/; 

SCALAR NODECAP NODE CAPACITY /1/; 

SCALAR REQTIM MAXIMUM TIME FOR THE OPERATION /25/; 


VARIABLES 
XCK,N,M,T) FLOW AMOUNT FROM N TO M FOR UNIT K AT TIME T 
Z TOTAL AMOUNT TIME TO TRANSPORT FROM SOURCE TO SINK NODE; 
BINARY VARIABLES X; 


EQUATIONS 
COST DEFINE OBJECTIVE FUNCTION 
MIDTIM(K) TIME TO STAY AT MIDPOINT FOR UNIT K 
mck, M,T) INTERMEDIATE NODE FOR UNIT K AT TIME T 
STARTCK,M,T) SOURCE NODE N FOR UNIT K AT TIME T 
TSINK(K,N) SINK NODE N FOR UNIT K 
CAP(M,T) NODE CAPACITY FROM NODE N TO NODE M AT TIME T 


CAPACITY(N,M,T) ARC CAPACITY FROM NODE N TO NODE M AT TIME T; 
START(K,M,T)$(S(K,M,T) EQ 1).. SUM(N,ACM,N,T)*X(K,M,N,T)) =E= 1; 
MID(K,M,T)$(S(K,M,T) EQ 0 AND ORD(M) NE CARD(M)).. 

SUM(N,AC(M,N,T)*X(K,M,N,T))SC(ORD(T) NE REQTIM) 
-SUM(N,A(N,M,T-1)*X(K,N,M,T-1)) =E= 0 ; 


MoENKCK ,N)SCORDCN) EQ CARDCN)):. <SUM((M,1),ACM,N,T-1)*AX(K,M,N,T<1)) 


=E= -1; 


CAPACITY(N,M,T)$C(ORD(M) NE ORD(N)).. SUM(K,X(K,N,M,T)*ACN,M,T)) =L= 
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ARCCAP; 


CAP(M,T)$(CARD(M) NE ORD(M)).. SUM((K,N),ACN,M,T)*X(K,N,M,T)) =L= 
NODECAP; 


MIDTIM(K).. SUM((N,M,T)S(MIDICK,N,M) EQ 1),X(K,N,M,T)*AC(N,M,T)) 
=G= MINTIM; 


COST. . SUMCCK,N,M,T)SCORD(CM) EQ CARD(M) ) 
,ORDCE) “ACK NM, U=1 ACN MN, eb — oe 


MODEL Rill Abi; 

OPTION LIMROW =0, LIMCOL = 0 , SOLPRINT = OFF ,ITERLIM = 10000; 
SOLVE RT1 USING MIP MINIMIZING 2; 

DISPLAY X.L,Z. L; 
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APPENDIX B. GAMS PROGRAM FOR MODEL 2 


* LP PROGRAM FOR FLEXIBLE COLUMN LENGTH 

* maj Dong Keun, Lee 

* Model : Multicommodity dynamic transportation network to minimize 
ve the transportation time 

/M1*M16/ 

iis 

JRIAKS/; 


et ot NODE 
Pettik PERIODS 
K UNIT troops 
ALIAS(M,N); 
ALIAS(T,S); 
TABLE A(M,N,T,S) ROAD NETWORKBETWEEN NODE M AND N AT TIME T TO S 
ae leo ol Go Onl On 17 eh olor Jn oa Le 


ale 
M2: 
ulby 
M4. 
tS. 
iG: 
yee 
ng: 


yee 
M6. 
My. 


M1 
(M2,M4,M5) 
(M3,M4) 
(M5 ,N4) 
(M7 ,M5,M6) 
M8 


Cil6,M7 >) 
Choe M13) 
a 


(M7 ,MS ,M6) 
M8 
(M6 ,M7) 


i 
i 


TG 


1 


1 
i 


i 


a 


1 


1 


] 


if 


Moar? M13) 
a 


M6. M8 
M7. (M6,M7) 
M10. (M12,M10) 
M13.M15 
M9. (M9 ,M13) 

+ 


MiG aC M10} 
M12. (M12 ,M14) 
ule yy Bs 

M14. (M14 ,M15) 


— bb 


Til Tis tia ie hile 


1 


b+ + bs be 

b+ b+ b+ 

+ bb 
-— 


1 
1 
1 


l 


EL 28 


1 
1 


+ iho tert isalo .4. 16-05. 17 Shon 7. To er PLO i 


(ey 1 1 1 : il 1 N i 1 
MYO. M11 1 1 1 
M729 f 1 1 i 1 1 I 
+ LOS Pa2e tr ivi To Tia ioe Gls 
MIOeMi11 1 1 1 1 
M7. M9 1 il 
+ BA BF ao WS raed sre WC 88) Ss eon BLO es G1 Oe Wa 1D sane U6 a 
M8. M10 i 1 1 1 a i 1 
+ Tiss liy ois. rs 
Meer 12 1 i 
use to Lo. io) TIO. Diss lt. tis Tie Dloetls. bly ie. TAs 
fo. M9 1 il 1 i 1 1 
+ foeno 14h b> PhO MG yi ilely 1 12. 1es0 bool oe I ke 
M6. M6 1 i a t 1 1 t 
M8.M14 1 
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ay EE Ey (a 1 


Ws yal y ike: 1 i 1 1 
+ T10.T15 THe T16 Ti2zs0i7 tis ieee hl 
M8.M14 1 1 1 a 1 
ES ty 1 
Mi M16 1 1 1 
Mises if 1 1 1 Il 
+ T15. 720 T16.T21 717.922 Tleeiizea7 19° 124 
M1I5.M16 1 1 IL iL 1 1 
E Tiged24 Tia. 125 
M11.M16 J! 1 ; 
PARAMETERS 


SR(K,M,T,S) SOURCE NODE FOR UNIT K AT NODE M FROM T1 TO T2 
{eM Tia 
Keerab ce Ll ee) mel 
Koos. le ay, 

NARCCAP(M,N) ARC CAPACITYBETWEEN M AND N 


/M1. M7 60 
M1.M1 60 
M2.(M2,M4) 30 
M2. M5 30 


MoacMsr ia). ~30 
M4.(M4,M5) 30 
Maas My Mowe = 30 
Meat M8) 30 


M6. M6 150 
fe2MS 120 
als (OUR), 60 
gree Ws 30 
M7. M6 30 
M8.M8 30 
M8.M14 150 
ye) WOH) 90 
biome M11) 30 
Mot 0 
ig by Gales 30 


MLO... ( M11) 60 
MiO 7 COMO; Mi 2 30 
Mola = CMeiaa) 150 


M11.M12 90 
M11. M16 330 
M12. M12 30 
M12. M14 30 
M13.M13 150 
M13. M15 30 
M14.(M14,M15) 30 
M15.M15 30 
M15. M16 150/ 


MP(N) | MIDPOINT NODE 
/M6 1 
M11 1 
M13 1/ 


S1(K,M,T) SOURCE NODE FOR UNIT K AT TIME Tl 

/K1.M1.T1 1 
KOuM2 ee |. 
K3.M3.T1l 1/ 

NODEI(M) M IS SOURCE SUPPLY POINT SINK NODE 
/M1 1 
M2 1 
M31 
M6 1 

M11 1 

M13 1 

M16 1 / 


MIDPT(K,M) MIDPOINT NODE FOR UNIT K 

7A ngs eo a E 
helt 7 
Rov MOe s/s 

SCALAR TROOP NUMBER OF VEHICLE PER TROOP /30/; 

SCALAR NNODECAP CAPACITY OF NODE 7 30/; 

VARIABLES 

XCK,M,N,T,5) FLOW AMOUNT FROM M TO N TIME BETWEEN T TO S FOR UNIT K 


Z TOTAL AMOUNT FLOW * TIME : 
POSITIVE VARIABLE x; 
EQUATIONS 
COST OBJECTIVE FUNCTION 
SlanickK,N,1,S) SOURCE NODE FOR UNIT K AT NODE N TIME T TO S 
INTERMID(K,N,T) INTERMID NODE FOR UNIT K AT NODE N TIME T 


TSINK(CK,N) SINK NODE N FOR UNIT K 

MIDPOINT(K,N) MIDPOINT NODE N FOR UNIT K 

STAY(K,N,T) TIME TO STAY AY MIDPOINT NODE N FOR UNIT KAT TIME T 
ARCCAP(M,N,T,S ) ARC CAPACITY M TO N TIME T TO §S 

NODECAP(M,T) NODE CAPACITY M AT TIME T; 


StARTCK,N, 1,5 )S(SROK,N,T,S ) EQ 1).. SUM(M,ACN,M,1T,S)*X(K,N,M,T,S)) =E= 
TROOP; 
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INTERMID(K,N,T)$(S1(K,N,T) NE 1 AND ORD(N) NE CARD(N)).. 
See.) ACM Nes, 0) *XCh,MN,o,1)) + SUMCCN,S), ACN MT, S)*ACK NMS T,S 
)) =E= 0 ; 


TSINK(K,N)S$(ORD(N) EQ CARD(N)).. -SUM((M,S,T),A(M,N,S,T)*X(K,M,N,S,T)) 
=E= -TROOP; 


ARCCAP(M,N,T,S)SCMP(N) NE 1).. SUMCK,ACM,N,T,S)*X(K,M,N,T,S)) =L= 
NARCCAP(M,N); 


MIDPOINT(K,N)S(MIDPT(K,N) EQ 1).. SUM((M,T,S)S(ORD(M) NE ORD(N)) 
vA(M,N,T,S)*X(K,M,N,T,S) )=G= TROOP; 


NODECAP(M,T)SCNODEICM) EQ 0).. SUMCCK,S,N), ACN,M,S,7)*X(K,N,M,5S,7)) 
=L= NNODECAP; 


STAY(K,N,T)S(MIDPT(K,N) EQ 1).. SUM((M,S)S(ORD(M) NE ORD(N)),A(M,N,S,T) 


LS) 


*XCK,M,N,S,T)) =L= SUM(S,A(N,N,T,S)*X(K,N,N,1T,S)); 


COST.. SUM((K,N,M,T,S)$(ORD(M) EQ CARD(M)),ORD(S)*X(K,N,M,T,S)*ACN,M, 
eso) Sa ee 


MODEL CH3 /ALL/; 

OPTION ITERLIM = 10000 ,LIMROW = 0 ,LIMCOL = 0 , SOLPRINT = OFF ; 
SOLVE CH3 USING LP MINIMIZING Z ; 

DISPLAY X%.UgaG; 
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